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SYSTEMS AND METHODS FOR LOOP LENGTH AND BRIDGED TAP LENGTH 
DETERMINATION OF A TRANSMISSION LINE 

Field of the Invention 

This invention relates to determination of transmission line characteristics. In 
particular, this invention relates to systems and methods for determining loop lengths and 
bridged tap lengths of a transmission line. 

Background of the Invention 

The collection and exchange of diagnostic and test information between transceivers 
in a telecommunications environment is an important part of a telecommunications, such as 
an ADSL, deployment. In cases where the transceiver connection is not performing as 
expected, for example, where the data rate is low, where there are many bit errors, or the like, 
it is necessary to collect diagnostic and test information from the remote transceiver. This is 
performed by dispatching a technician to the remote site, e.g., a truck roll, which is time 
consuming and expensive. 

In DSL technology, communications over a local subscriber loop between a central 
office and a subscriber premises is accomplished by modulating the data to be transmitted 
onto a multiplicity of discrete frequency carriers which are summed together and then 
transmitted over the subscriber loop. Individually, the carriers form discrete, non- 
overlapping communication subchannels of limited bandwidth. Collectively, the carriers 
form what is effectively a broadband communications channel. At the receiver end, the 
carriers are demodulated and the data recovered. 

DSL systems experience disturbances from other data services on adjacent phone 
lines, such as, for example, ADSL, HDSL, ISDN, Tl, or the like. These disturbances may 
commence after the subject ADSL service is already initiated and, since DSL for internet 
access is envisioned as an always-on service, the effect of these disturbances must be 
ameliorated by the subject ADSL transceiver. 

SUMMARY OF THE INVENTION 

Identifying, measuring and characterizing the condition of a transmission line is a key 
element of an ADSL deployment. In cases when the transceiver connection is not performing 
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as expected, for example, the data rate is low, there are many bit errors, a data link is not 
possible, or the like, it is important to be able to identify the loop length and the existence, 
location and length of any bridged taps without having to send a technician to the remote 
modem site to run diagnostic tests. 
5 This invention describes a system and method for estimating the loop length, the 

number of bridged taps and length of the bridged taps on a transmission line from readily 
available modem data. The loop length, the number of bridge taps and the length of the 
bridged taps can be estimated by comparing a measured frequency domain channel impulse 
response of the transmission line to a model of a transmission line that is composed of 

10 multiple sections and multiple bridge taps. The diagnostic and test information describing the 
condition of the line can then be exchanged, for example, by two transceivers during a 
diagnostic link mode, such as that described in copending Attorney Docket No. 
081513.00004, filed herewith and incorporated herein by reference in its entirety. 

These and other features and advantages of this invention are described in or are 

15 apparent from the following detailed description of the embodiments. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The embodiments of the invention will be described in detail, with reference to the 
following figures wherein: 
20 Fig. 1 illustrates an exemplary multiple section loop with multiple bridged taps; 

Fig. 2 illustrates a graph of the measured received reverb signal and the theoretical 
model for downstream data; 

Fig. 3 illustrates a graph of the measured received reverb signal and the theoretical 
model for upstream data; 
25 Fig. 4 is a functional block diagram illustrating an exemplary loop length and bridged 

tap length estimation system according to this invention; 

Fig. 5 is a flowchart outlining an exemplary general method for determining loop 
length and bridged tap lengths according to this invention; 

Fig. 6 is a flowchart outlining an exemplary method for estimating the loop length and 
30 bridged tap length in the upstream direction according to this invention; and 

Fig. 7 is a flowchart outlining an exemplary method for estimating the loop length and 
bridged tap length in the downstream direction according to this invention. 
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DETAILED DESCRIPTION OF THE INVENTION 

The exemplary embodiments of this invention will be described in relation to the 
application of the invention to an ADSL transceiver environment. However, it should be 
appreciated that in general the systems and methods of this invention will work equally well 
5 for any multiple section loop with one or more bridged taps. 

For example, during the ADSL modem initialization, the frequency domain channel 
impulse response of the subscriber loop is measured at a set of discrete frequency values. 
The measured frequency values are designated as H m (/j), and / = i A/ for i = 0,1, k-1, 
where A/is the frequency spacing between adjacent samples. 
10 Fig. 1 illustrates an exemplary model of a loop with N sections and M bridged taps. 

The frequency domain model for the channel impulse response for the loop in Fig. 1 can be 
written as H(x/), where /is the frequency and the vector x contains the lengths (d.) of the N 
sections of the loop and the lengths (bj)of the M bridged taps: 



15 x = [di d2....,d N ,b ! ,b 2 ;- .,bM,]- 

Assuming that the number of sections of the multiple section subscriber loop, N, and the 
number of bridged taps, M, are known, an estimate of the optimal parameter vector x that 
best approximates the measured channel impulse response H m (//) can be determined given the 
20 model The optimal parameter vector set x* can be estimated by minimizing the norm 

of the difference between the measured and the model frequency response, at the discrete 
frequency values / = z A/ for i - 0,1, k-1. This minimization can be performed using the 
expression: 

x* = min2|H m (^)-H(x,/:)f 2 . 

X i=0 



25 If the number of the bridged taps on the loop is not known, by adopting a large number of 
bridged taps in the model frequency response, and assuming that the minimization will 
converge to a solution with the correct number of bridged taps with non-zero length, the 
remaining bridge taps will have length zero. 

The frequency domain model H(x/) can also incorporate the effect of, for example, 

30 an imperfectly matched transmission line, by including the effects of the load and source 
impedances. 

3 
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More particularly, the loop characterization algorithms employ a model based 
approach to estimate the length of the loop and the lengths of up to two bridged taps. A 
channel characterization algorithm compares the measured channel impulse response to the 
channel impulse response of a loop model consisting of a single-gauge wire and containing 
5 up to two bridged taps. However, it is to be appreciated that the basic model can be extended 
to include multiple gauge wires and multiple bridged taps. The loop length and the bridged 
tap lengths are the parameters of the theoretical channel impulse response. The system varies 
the parameters of the theoretical model and evaluates the difference between the measured 
channel impulse response and the theoretical channel impulse response. The loop 
1 0 length/bridged tap lengths that minimize the error function are then declared as the estimated 
values. The presence of a bridged tap is declared if the bridged tap length is greater than a 
predetermined length, such as one hundred feet. This threshold for bridged tap detection was 
set experimentally. It was determined that for most loops there is a chance that a phantom 
bridged tap with a small length will be detected because of modeling inaccuracies and noise 
1 5 in the measurement system. Since the lengths of these phantom bridged taps were almost 
always below 100 ft, the exemplary threshold was set to 100 ft. However, in general the 
threshold can be altered depending on the particular operational environment and the 
complexity of the model. 

There are two separate algorithms which perform loop characterization for 
20 downstream (DS) and upstream (US) data. For example, during modem initialization, data 
collection software collects the reverb signal by averaging K consecutive frames where K > 
64. However, it is to be appreciated that as more averaging is performed, the less noisier the 
measurement will be. However, since there is a prescribed number of frames in the standard 
modem training where the reverb signal is transmitted, the exemplary number of averages 
25 was set at 64. The received reverb signal obtained in this way is an estimate of the impulse 
response of the entire channel including the front-end responses of the transmitting and 
receiving modems. The frequency domain received reverb signal is obtained in accordance 
with: 

30 where f is a dummy variable denoting frequency and rx(n), for n = 1 , . . . , N, are the samples 
of the time-domain received reverb signal within a frame, with N being the number of 
samples contained in a single frame. Equation 1 may contain a slight abuse of notation 
because in reality the frequency variable f is not continuous but rather discrete and for this 
4 
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reason the channel impulse response is available at a set of discrete frequencies, called 
tones, which are multiples of Af- 4312.5Hz: 



fi = iAf,i=\,...m. (2) 

The reverb signal is transmitted over a portion of the entire ADSL spectrum. For example, 
5 the reverb signal is available at 224 (96 in GXite) tones from fa = 32 Af to fas = 255A/in the 
downstream channel and at 26 tones from/6 = 6Af to fa = 31 A/*in the upstream channel. The 
downstream reverb signal is collected at the customer-premises equipment (CPE) and 
upstream reverb signal is collected at central office (CO). While there is no difference in the 
data collection process for the upstream or the downstream reverb signal, the characteristics 

10 of these two data sets are quite different. Specifically, the downstream reverb data contains 
significantly more information. Furthermore, there are more samples of the frequency 
domain reverb signal available in the downstream direction and these samples cover an 
extended range in the frequency domain where the effects of bridged taps on impulse 
response can be easily detected. However, there is one crucial difference between the 

15 upstream and the downstream data sets which complicates using the same interpretation 

algorithm for both. In the downstream channel, the matching of the front-end impedance to 
the loop impedance tends to be better than in the upstream channel. This makes it possible to 
use a simplified channel model for the downstream channel. Unfortunately, the impedance 
matching in the upstream channel is generally not as good as in the downstream channel and 

20 a more complicated channel impulse response should be used. 

Due to these complications in channel modeling, and the lack of sufficient data 
samples, the basic upstream channel characterization algorithm is limited in terms of 
estimation accuracy and the number of bridged taps that can be detected. However, by 
extending the channel model to include multiple sections of varying gauges and/or more than 

25 two bridged taps, the presence of more than two bridged taps can be detected and more 
accurate results for the lengths of individual sections of the loop determined if there is a 
change of wire gauge along the loop. The only trade off is that as the number of model 
parameters increase, the computational effort needed to estimate the parameters will increase 
as well. 

30 The following describes the theoretical details leading to the derivation of the 

frequency domain channel impulse response of the model and explains the channel 
characterization for both the downstream and the upstream data in detail. Both the 
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downstream and the upstream interpretation algorithms employ the same least squares 
minimization concept where the square of the error norm between the actual and the 
theoretical channel impulse responses is niinimized, but differ in the theoretical channel 
impulse response used. 

For the loop characterization for downstream data, an exemplary two-wire loop is 
characterized by its characteristic impedance: 



+ jaL 
jaC 



And its propagation constant: 



r(f) = + j<t*%P + J.* 00 ) 
1 0 where co = 2nf is the radian frequency and R (resistance), L (inductance), G (admittance) and 
C (capacitance) are the frequency dependent constants of the loop and vary with wire gauge. 
For a perfectly terminated loop, or a very long loop, with length d, and two bridged taps of 
lengths bi and b 2 , the transfer function of the loop H(d, b], b 2 ,f), is given by: 

e -dr(J) 

H(rf,6 1 >,/) = [2 + tanh(M][2 + tanh(62; , )] (3) 

15 In logarithmic scale: 

log|H(d,bi,b 2l #| = log(2)-dy(/) - log[2+tanh(b,y)] - log[2+tanh(b 2 y)]. (4) 

Note the linear dependence of the loop loss to the length of the cable. The actual 
transfer function of the loop can be measured during modem initialization. . Then the 
measured transfer function of the loop is matched with that of a loop of length d with two 
20 bridged taps as given in Eq. 3 . In other words, determining d, b u and b 2 minimizes the 
following least squares error criterion: 

tmnt\H(d,b h b2,f!)-Rx(M ( 5 ) 
d ' b - b /=* 

25 where Rxffi) is the received reverb signal sampled at^ =i/and z',- are U are the first and the 
last tones Rx{fl). 
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An example of the operation of the algorithm for an exemplary loop is illustrated in 
Fig. 2. Displayed are the measured received reverb signal Rx{f) and the theoretical model H 
(d, b\,b 2 j) which were obtained by finding the model parameters d,b\,b 2 that best match the 
data. Specifically, the observed (dashed line) received reverb signal Rx(f) is plotted against 
5 the theoretical channel model (solid line) HidMMJ) as functions of frequency for an 

exemplary 6000 ft loop with an exemplary single 1300 ft bridged tap. The exemplary loop 
consisted of a 26 awg. 6000 ft wire with a 26 awg. 1300 ft bridged tap close to the CPE. The 
model parameters best matching the observed data were found to be d=6000 ft, b\ = 1300 ft 
and b 2 = 0 ft. 

10 It follows from Eq. 5 that the interpretation algorithm basically does a search over the 

variables J, b\ and b 2 and finds the ones minimizing the cost function given below: 

E{d,bxM = £| H(d,biMfi) - Rx(JI) g. (6) 

i=!f 

15 Since the cost function E(dM,b 2 ) is a nonlinear of d, b x and b 2 , the function contains many 
local minima. Therefore, many well known optimization algorithms such as Gauss-Newton 
should not be used since these algorithms are unable to cope with multiple local minima and 
they converge to a local minimum of the cost function. In this exemplary embodiment the 
global minimum of E{d,b\M) is desired. For this reason, a brute-force global minimization 

20 algorithm is used where the cost function is sampled at the points (d p ,tf ,b r 2 ) , (f = pAD, 
b\ =qAbi and b r 2 = rAbz with;? = l,...,P, q = 1.....Q and r= 1,...,R. Next the parameters 
{d p ,tf ,b r 2 ) which result in the minimum cost among the sampled values are chosen. This 
requires evaluating the cost function at P x Q x R locations. 

In order to be able to determine the theoretical transfer function of the loop, 

25 H(d, b\, bi,f), the frequency dependent propagation constant y(f) for a number of wires of 

different gauges needs to be stored. In an exemplary embodiment 24 awg. and 26 awg. wires 
are used which require 4 x N locations to store the real and the imaginary parts of y(J) for N 
ADSL tones. Additionally, the analog front end (AFE) compensation curves need be stored 
which occupy N locations in memory. Depending on where the algorithm is implemented, 

3 0 the loop transfer function can be determined directly from Eq. 4, for example, if the algorithm 
were implemented on a personal computer or workstation, or it may be necessary to store the 
log[2+tanh(6iY)] terms in regular intervals as required by the sampling procedure for 
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(d p , b?,b r 2 ) . For example, it is possible to pre-compute and store the log[2 + tanh(6,Y)], i = 
1,2, from b\ = 100 ft to b\ - 2000 ft to in 100 ft intervals. Assuming low processor power, 
the log[2 + tanh(Z?iY)] terms can be predetermined and stored which take about 20 x N 
locations for the real part only. Therefore, in this exemplary embodiment, the total memory 
5 is about (20+4+1+3) x N = 28 x N where 2 x 256 locations are needed to store intermediate 
variables determined during the execution of the algorithm. 

Although it will not be shown here, it is possible to simplify the computation of the 
cost function E(d,bj,b 2> ) so that only 12 multiplications and 15 additions are needed. This 
means that the total computational complexity of the algorithm is about PxQxRx(/i 

10 multiplications + 15 additions) plus some additional start-up computations which are 
negligible compared to the above figure. 

Unlike the downstream interpretation case, for upstream interpretation it is more 
accurate to assume that the line is not perfectly terminated. Specifically, the impedance 
mismatch at the transmitter-line connection at the CPE modem and the impedance mismatch 

15 at the receiver-line connection at the CO modem become important factors that should be 
taken into account. While the basic idea behind the channel characterization algorithm for 
the upstream data remains the same, and involves matching a theoretical channel transfer 
function to the actual measured transfer function, the computation of the theoretical channel 
transfer function becomes much more involved. As with the downstream interpretation case, 

20 the channel transfer function is again measured by averaging K frames of the received reverb 
signal as given by Eq. 1 . 

The theoretical model for the channel transfer function in the upstream case can be 
described in two steps. The first step consists of writing the equations for the current and the 
voltage at the source (CPE), I s , V s , in terms of current and voltage at the load (CO), I L , V u 

25 through the application of ABCD matrices: 



where A', B, F s and F L are 2 x 2 matrices whose elements are arrays of N elements. Here, A' 
is a matrix representing the frequency domain response of the z'th section of the loop, B is the 



representing the frequency domain response of the analog front end (AFE) hardware of the 




(7) 



30 matrix representing the response of the bridged tap and F s and F L are the matrices 
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modem circuitry for TX (source) and RX (load) paths. From Eq. 7 the transfer function of 
the channel can be derived and is given by: 

H(d u d 2 ,b,f) = ^, (8) 

5 

where d\ is the length of the section before a bridged tap and d 2 is the length of the section 
after the bridged tap. Note that the CO interpretation algorithm uses a two-section, single 
bridged tap model. This is because of the limited number of frequency bins,/i=i'4f, from 
tone i=6 to i=32, at which the transfer function is available. 
1 0 Entries of the above matrices are given as follows: 

4', = 4a = coshOc/,.) 

4=z 0 sinhorf,), 4=4 2 z 0 - 2 

Entries of matrix B: 

15 B n =0,B 2l =Z] l (b) 
Where ZJ 1 =tanh(b^/Z 0 , and finally: 

F n =F 22 -=1, F£ = 0, F^-ZJ 1 
The estimation algorithm minimizes the difference between the measured and the actual 
20 transfer functions: 

r^|//(4,rf 2 A/)-W)E. (9) 

An example of the operation of the upstream loop length and bridged tap length 
25 estimation algorithm is illustrated in Figure 3. Here the measured received reverb signal 
Rx(f) and the theoretical model H(d,b\.b 2 ,f), which was obtained by finding the model 
parameters d,b\b 2 that best match the data, are displayed. The exemplary loop consisted of 
26 awg. 7700 ft wire with a 26 awg. 600 ft bridged tap 5900 ft away from CO. The model 
parameters best matching the observed data were found to be di=7900 ft, d 2 = 0 ft and b = 
30 500 ft. Note that although the d x and d 2 parameters found by the algorithm are different than 
9 
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their actual values, the actual values are d } = 5900 ft and d 2 ~ 1800 ft, the sum of d\ + d 2 is 
within 200 ft of the actual loop length. This example illustrates that even though the loop 
length is fairly accurate the location of the bridged tap is difficult to reliably estimate. 

From the expressions leading to the theoretical channel transfer function, H(di,d 2 ,b,f), 
5 it is clear that for the exemplary computation of the theoretical channel response Z s , Z L , Z 0 
and y, for 24 awg. and 26 awg., need be stored, and that characterizing the bridged tap 
is dependent on the bridged tap length. Assuming an exemplary resolution of 100 ft in 
bridged tap length and a maximum exemplary detectable bridged tap length of 2000 ft, there 
are 20 different Z/6i) arrays. Finally, the sinh(.) and cosh(,) elements of the matrices A] and 

10 A2 are stored. Then, assuming a 500 ft resolution in loop length and a maximum measurable 
loop length of 20,000 ft, there should be 80 x 46 locations for storing entries of A,-. In total 
for storing these variables there should be 108 x 46 memory locations, including storage for 
Rx(j) and H(d x ,d 2 ,b,f), and another 10 x 46 locations are needed for storing intermediate 
variables during the execution of the algorithm, giving a total of approximately 1 1 8 x 46 

1 5 memory locations for this exemplary embodiment. 

Figure 3 illustrates the observed (dashed line) received reverb signal Rx(f) plotted 
against the theoretical channel model (solid line) H (d u d 2 ,f) as functions of frequency for an 
exemplary 7700 ft loop with a single 600 ft bridged tap. 

During the search process, P values for d\ , Q values for b and R values for d 2 are 

20 selected and the cost function for each combination of d\,di,b determined. Thus, to determine 
the channel impulse response there are 4 x (8 x 23 complex multiplications + 4 x 26 complex 
additions). Therefore the total computational cost in this exemplary embodiment is P x Q x R 
x (32 x 26 complex multiplications + 4 x 26 complex additions). 

Fig. 4 illustrates an exemplary loop length and bridged tap length estimation system 

25 according to an embodiment of this invention for downstream data. In particular, the loop 
length and bridged tap length estimation system 100 comprises a downstream loop length and 
bridged tap length determination device 200, an upstream loop length and bridged tap length 
determination device 300, a central office modem 20 and a consumer-premises modem 30, 
connected by link 10, such as a twisted pair. The a downstream loop length and bridged tap 

30 length determination device 200 comprises a controller 210, an I/O interface 220, a storage 

device 230, a reverb signal determination device 240, a loop length output device 250 and a 

bridged tap output device 260, connected by link 5. The upstream loop length and bridged tap 

length determination device 300 comprises a controller 310, an I/O interface 320, a storage 

device 330, a reverb signal determination device 340, an impedance determination device 
10 
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350, a modem identification device 360, a loop length output device 370 and a bridged tap 
output device 380, connected by link 5. 

While the exemplary embodiment illustrated in Fig. 4 shows the components of the 
loop length and the bridged tap length estimation system and associated components 
5 collocated, it is to be appreciated that the various components of the loop length and the 
bridged tap length estimation system 100 can be located at distant portions of a distributed 
network, such as a local area network, a wide area network, an intranet and/or the Internet, or 
within a dedicated loop length and bridged tap length estimation system. Thus, it should be 
appreciated that the components of the loop length and bridged tap length estimation system 

10 100 can be combined into one device or collocated on a particular node of a distributed 
network. As will be appreciated from the following description, and for reasons of 
computational efficiency, the components of the loop length and the bridged tap length 
estimation system 100 can be arranged at any location, such as in a general purpose computer 
or within a distributed network without affecting the operation of the system. 

15 Furthermore, the links 5 can be a wired or a wireless link or any other known or later 

developed element(s) that is capable of supplying electronic data to and from the connected 
elements. 

In operation, for determination of the loop length and the bridged tap length in the 
downstream direction, the controller 210, in cooperation with the I/O interface 220 triggers 

20 initialization of the modem 20. The reverb signal determination device 240, in cooperation 
with the modem 20, the controller 2 1 0 and the I/O interface 220 determines a transfer 
function by averaging K consecutive frames of a reverb signal. The loop length, a first 
bridged tap length and a second bridged tap length are input from an input device (not shown) 
such as a computer, a laptop, a terminal, a transmission line testing device, or the like, or 

25 retrieved from the storage device 230. 

The controller 210, in cooperation with the storage device 230, then determines the 
frequency domain propagation function for a specified wire gauge, and the frequency domain 
loop model. The calibrated and compensated reverb signals in the frequency domain are 
stored in the storage device 230 and the reference wire gauge input or retrieved from the 

30 storage device 230. 

The controller 210, in cooperation with the storage device 230 determines the number 

of elements in the Rx function and the difference between the actual and the measured 

transfer function. The loop length output device, in cooperation with the I/O interface then 

outputs the estimated loop length to, for example, a computer, a laptop, a terminal, a 
11 
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transmission line testing device, or the like. Additionally, the bridged tap output device 
outputs the estimated bridged tap length to, for example, a computer, a laptop, a terminal, a 
transmission line testing device, or the like. 

In operation, for detennination of the loop length and bridged tap length in the 
5 upstream direction, the controller 3 1 0, in cooperation with the I/O interface 320 triggers 
initialization of the modem 30. The reverb signal determination device 340, in cooperation 
with the modem 30, the controller 310 and the I/O interface 320 determines a transfer 
function by averaging K consecutive frames of a reverb signal. 

Next, the controller 310, in cooperation with the storage device 230, determines the 
10 frequency domain propagation function for a specified wire gauge, where the specified wire 
gauge is input or retrieved from the storage device 330. 

The controller 310, in cooperation with the storage device 330 and the impedance 
detennination device 350, determines the frequency domain impedance of the specified wire 
gauge. Then, the controller 310, in cooperation with the storage device 330 and the 
15 impedance determination device 350, determines the transmit impedance of the CPE modem 
and the receive impedance of the CO modem. 

The controller 310, in cooperation with the storage device 330, determines the matrix 
representing the frequency domain responses of the i ,h section of the loop, the matrix 
representing the response of the bridged tap, and the F s matrix representing the AFE circuitry 
20 for the source (TX) and load (RX) paths and stores them in the storage device 330, and 
estimates the transfer function H. The calibrated and the compensated reverb signal in the 
frequency domain and the reference gauge of the wire are input or retrieved from the storage 
device 330. 

The modem identification determining device 360 then determines the identification 
25 of the CO modem collecting the upstream reverb signal, and the identification of the CPE 
modem transmitting the upstream reverb signal. Knowing the number of elements in the Rx 
function, the controller 310 minimizes the difference between the actual and measured 
transfer functions, and outputs, with the cooperation of the loop length output device 370 and 
the bridged tap output device 380, the estimated loop length and the estimated bridged tap 
30 length, respectively. 

Fig. 5 illustrates an exemplary method of deterrmning a loop length and bridged tap 
lengths. In particular, control begins in step SI 00 and continues to step SI 10. In step SI 10, 
the channel impulse response is estimated based on a measured reverb signal. Next, in step 
S120, the theoretical channel impulse response of a loop model is determined using a loop 



BNSDOCID: <WO_ 



_0152439A1_I_> 



WO 01/52439 PCT/US01/00420 
length and the bridged tap lengths. Then, in step S130, the loop length and the bridged tap 
lengths of the model are varied. Control then continues to step S 140. 

In step SI 40, the difference between the measured channel impulse response and the 
theoretical channel impulse is monitored. Next, in step SI 50, the estimated values of the loop 
5 length and bridged tap length are declared based on the loop lengths and bridged tap lengths 
that minimize the error function between the measured channel impulse response and the 
theoretical channel impulse response. Control then continues to step S 160 where the control 
sequence ends. 

Fig. 6 illustrates an exemplary method of determining the loop length and the bridged 

10 tap length for downstream data. In particular, control begins in step S200 and continues to 
step S210. In step S210, a modem is initialized. Next, in step S220, a transfer function is 
determined by averaging K consecutive frames of the reverb signal. Then, in step S230, the 
loop length is input. Control then continues to step S240. 

In step S240, a first bridged tap length is input. Next, in step S250, a second bridged 

15 tap length is input. Next, in step S260, the frequency domain propagation function is 
determined for a specified wire gauge. Control then continues to step S270. 

In step S270, the frequency domain loop model is determined. Next, in step S280, the 
calibrated and compensated reverb signals in the frequency domain are input. Then, in step 
S290, the reference wire gauge is input. Control then continues to step S300. 

20 In step S300, the number of elements in the Rx function are input. Next, in step S3 1 0, 

the difference between the actual and the measured transfer function are determined. Then, 
in step S320, the estimated loop length is determined. Control then continues to step S330. 

In step S330, the estimated bridged tap length is determined. Control then continues 
to step S340 where the control sequence ends. 

25 Fig. 7 illustrates an exemplary method of determining the loop length and bridged tap 

length for upstream data. In particular, control begins in step S500 and continues to step 
S5 10. In step S5 10, the modem is initialized. Next in step S520, the transfer function is 
determined by averaging K consecutive frames of the reverb signal. Then, in step S530, the 
frequency domain propagation function for the wire gauge in use is determined. Control then 

30 continues to step S540. 

In step S540, the frequency domain impedance of the wire gauge is determined. 
Next, in step S550, the transmit impedance of the CPE modem is determined. Then, in step 
S560, the receive impedance of the CO modem is determined. Control then continues to step 
S570. 

13 
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In step S570, the. matrix representing the frequency domain responses of the i th section 
of the loop are determined. Next, in step S580, the matrix representing the response of the 
bridged tap is determined. Then, in step S590, the F s matrix representing the AFE circuitry 
for the source (TX) and load (RX) paths are determined. Control then continues to step S600. 
5 In step S600, the transfer function H is estimated. Next, in step S610, the calibrated 

and the compensated reverb signal in the frequency domain are input. Then, in step S620, the 
reference gauge of the wire is input. Control then continues to step S630. 

In step S630, the identification of the CO modem collecting the upstream reverb 
signal is input. Next, in step S640, the identification of the CPE modem transmitting the 

10 upstream reverb is input. Then, in step S650, the number of elements in the Rx function are 
input. Control then continues to step S660. 

In step S660, the difference between the actual and measured transfer functions are 
minimized. Next, in step S670, the estimated loop length is determined. Then, in step S680, 
the estimated bridged tap length is determined. Control then continues to step S690 where 

15 the control sequence ends. 

As illustrated in FIG. 4, the loop length and bridged tap length estimation system can 
be implemented either on a single program general purpose computer, or a separate program 
general purpose computer. However, the loop length and bridged tap length estimation 
system can also be implemented on a special purpose computer, a programmed 

20 microprocessor or microcontroller and peripheral integrated circuit element, an ASIC or other 
integrated circuit, a digital signal processor, a hard wired electronic or logic circuit such as a 
discrete element circuit, a programmable logic device such as a PLD, PLA, FPGA, PAL, a 
modem, or the like. In general, any device capable of implementing a finite state machine 
that is in turn capable of implementing the flowcharts illustrated in FIG. 5-7 can be used to 

25 implement the loop length and bridged tap length estimation system according to this 
invention. 

Furthermore, the disclosed method may be readily implemented in software using 
object or object-oriented software development environments that provide portable source 
code that can be used on a variety of computer or workstation hardware platforms. 
30 Alternatively, the disclosed loop length and bridged tap length estimation system may be 
implemented partially or fully in hardware using standard logic circuits or VLSI design. 
Whether software or hardware is used to implement the systems in accordance with this 
invention is dependent on the speed and/or efficiency requirements of the system, the 
particular function, and the particular software or hardware systems or microprocessor or 
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microcomputer systems being utilized. The loop length and bridged tap length estimation 
systems and methods illustrated herein, however, can be readily implemented in hardware 
and/or software using any known or later-developed systems or structures, devices and/or 
software by those of ordinary skill in the applicable art from the functional description 
5 provided herein and a general basic knowledge of the computer arts. 

Moreover, the disclosed methods may be readily implemented as software executed 
on a programmed general purpose computer, a special purpose computer, a microprocessor, 
or the like. In these instances, the methods and systems of this invention can be implemented 
as a program embedded on a personal computer such as a Java® or CGI script, as a resource 

10 residing on a server or graphics workstation, as a routine embedded in a dedicated loop length 
and bridged tap length estimation system, a modern,^ dedicated loop length and/or bridged 
tap length estimation system, or the like. The loop length and bridged tap length estimation 
system can also be implemented by physically incorporating the system and method into a 
software and/or hardware system, such as the hardware and software systems of a dedicated 

15 loop length and bridged tap length estimation system or modem. 

It is, therefore, apparent that there has been provided, in accordance with the present 
invention, systems and methods for loop length and bridged tap length estimation. While this 
invention has been described in conjunction with a number of embodiments thereof, it is 
evident that many alternatives, modifications and variations would be or are apparent to those 

20 of ordinary skill in the applicable arts. Accordingly, it is intended to embrace all such 

alternatives, modifications, equivalents and variations that are within the spirit and scope of 
this invention 
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What is Claimed is: 

1 . A loop length and bridged tap length estimation system comprising: 

a channel impulse response device that determines a channel impulse response 
for a transmission line; 

5 a modeling device that determines a theoretical channel impulse response; and 

a loop length and bridged tap determination device that determines a loop 
length and at least one bridged tap length based on a comparison of the channel impulse 
response and the theoretical channel impulse response. 

2. The system of claim 1 , wherein the comparison is based on an error 
10 minimization function. 

3. The system of claim 2, where the error miriimization function is a least squares 
minimization. 

4. The system of claim 1, further comprising an impedance determination device 
that determines at least one of a transmit impedance of a CPE modem and a receive 

1 5 impedance of a CO modem. 

5. The system of claim 1, further comprising a modem identification detennining 
device that determines at least one of an identification of a CPE modem and an identification 
of a CO modem. 

6. The system of claim 1, wherein the theoretical channel impulse response 
20 comprises variables representing a load and a source impedance. 

7. The system of claim 1 , where the channel impulse response is estimated based 
on a measured reverb signal. 

8. The system of claim 1, wherein the transmission line comprises at least one of 
at least one wire gauge and at least one bridged tap. 

25 9. The system of claim 1, wherein the transmission line is used for multicarrier 

modulation communications. 

10. A method for determining a loop length and at least one bridged tap for a 
transmission line comprising: 

estimating a channel impulse response; 
30 determining a theoretical channel impulse response; 

monitoring a difference between the theoretical channel impulse response and 
the estimated channel impulse response; and 

determimng the loop length and at least one bridged tap based on an error 
minimization function. 

16 
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1 1 . The method of claim 1 0, further comprising determining a transfer function of 
the transmission line based on a reverb signal. 

12. The method of claim 1 1 , wherein the transfer function is determined based on 
averaging a predetermined number of consecutive frames of the reverb signal. 

5 13. The method of claim 1 0, further comprising determining a frequency domain 

propagation function for at least one wire gauge. 

14. The method of claim 10, further comprising determining at least one of a 
calibrated and a compensated reverb signal in the frequency domain. 

1 5 . The method of claim 1 0, further comprising determining the frequency 
1 0 domain impedance of at least one wire gauge. 

16. The method of claim 10, further comprising determining a transmit impedance 
of a CPE modem. 

17. The method of claim 1 0, further comprising determining a receive impedance 
of a CO modem. 

15 18. The method of claim 1 0, further comprising detennining a matrix representing 

a frequency domain response for a predetermined portion of the transmission line. 

19. The method of claim 10, further comprising determining a matrix representing 
a response of the at least one bridged tap. 

.20. The method of claim 10, further comprising determining a matrix representing 
20 an analog front end circuit for a source and a load path. 

21. The method of claim 1 0, further comprising determining an identification of at 
least one of a CO modem and a CPE modem. 

22. The method of claim 10, further comprising outputting an estimated loop 
length of the transmission line. 

25 23. The method of claim 10, further comprising outputting an estimated bridged 

tap length for the transmission line. 

24. An information storage media comprising information that determines a loop 
length and at least one bridged tap for a transmission line comprising: 
information that estimates a channel impulse response; 
30 information that determines a theoretical channel impulse response; 

information that monitors a difference between the theoretical channel impulse 
response and the estimated channel impulse response; and 

information that determines the loop length and at least one bridged tap based 
on an error minimization function. 
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25 . The information storage media of claim 24, further comprising information 
that determines a transfer function of the transmission line based on a reverb signal. 

26. The information storage media of claim 25, wherein the transfer function is 
determined based on averaging a predetermined number of consecutive frames of the reverb 

5 signal. 

27. The information storage media of claim 24, further comprising information 
that determines a frequency domain propagation function for at least one wire gauge. 

28. The information storage media of claim 24, further comprising information 
that determines at least one of a calibrated and a compensated reverb signal in the frequency 

1 0 domain. 

29. The information storage media of claim 24, further comprising information 
that determines the frequency domain impedance of at least one wire gauge. 

30. The information storage media of claim 24, further comprising information 
that determines a transmit impedance of a CPE modem. 

15 31. The information storage media of claim 24, further comprising information 

that determines a receive impedance of a CO modem. 

32. The information storage media of claim 24, further comprising information 
that determines a matrix representing a frequency domain response for a predetermined 
portion of the transmission line. 

20 32. The information storage media of claim 24, further comprising information 

that determines a matrix representing a response of the at least one bridged tap. 

33. The information storage media of claim 24, further comprising information 
that determines a matrix representing an analog front end circuit for a source and a load path. 

34. The information storage media of claim 24, further comprising information 
25 that determines an identification of at least one of a CO modem and a CPE modem. 

35. The information storage media of claim 24, further comprising information 
that outputs an estimated loop length of the transmission line. 

36. The information storage media of claim 24, further comprising information 
that outputs an estimated bridged tap length for the transmission line. 

30 
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